Angular correlations in the cosmic gamma-ray background from 
dark matter annihilation around intermediate-mass black holes 



00 

o 
o 

(N 

> 
O 

(N 



I 

o 



> 

On 



oo 

o 



X 



Marco Taoso,^'^ Shin'ichiro Ando,'^ Gianfranco Bertone,^ and Stefano Profumo'* 

^INFN, Sezione di Padova, via Marzolo 8, Padova, 35131, Italy 
^Institut d Astrophysique de Paris, UMR 7095-CNRS, 
Universite Pierre et Marie Curie, 98bis boulevard Arago, 75014 Paris, France 
^California Institute of Technology, Mail Code 130-33, Pasadena, California 91125, USA 
^ Santa Cruz Institute for Particle Physics and Department of Physics, 
University of California, Santa Cruz, California 95064, USA 

Dark matter (DM) annihilation could in principle contribute to the diffuse cosmic gamma-ray back- 
ground (CGB). While with standard assumptions for cosmological and particle physics parameters 
this contribution is expected to be rather small, a number of processes could boost it, including 
a larger-than-expected DM annihilation cross-section, or the occurance of DM substructures such 
as DM mini-spikes around intermediate-mass black holes. We show that angular correlations of 
the CGB provide a tool to disentangle the signal induced by DM annihilation in mini-spikes from 
a conventional astrophysical component. Treating blazars as a known background, we study the 
prospects for detecting DM annihilations with the Fermi Gamma-Ray Space Telescope for different 
choices of DM mass and annihilation channels. 



I. INTRODUCTION 

The identification of nonbaryonic dark matter (DM) 
is still an open problem. Weakly interacting massive 
particles (WIMPs) are among the best motivated DM 
candidates, due to their connections with several, inde- 
pendently formulated particle physics theories beyond 
the Standard Model, and also in view of their intrigu- 
ing phenomenology (see Refs. [E 0, H, Q for reviews). 
WIMPs are actively searched for with underground de- 
tectors, with searches for related signatures with the 
Large Hadron Colider, and, indirectly, through the de- 
tection of their annihilation products. 

In particular, within the latter detection technique, an- 
timatter and gamma-ray signals from WIMPs annihila- 
tion are actively being searched for with the PAMELA 
(Payload for Anti-Matter Exploration and Light-nuclei 
Astrophysics) satellite and with the Fermi Gamma-ray 
Space Telescope (formerly known as GLAST), both cur- 
rently taking data. PAMELA has actually already found 
an interesting feature in the positron ratio spectrum 
Ref. [SL made even more interesting by the recent ATIC 
data [6|. This feature could be explained in terms of 
the annihilation or decay of DM particles, although the 
gamma-ray flux implied by these models severely con- 
strain this interpretation (see Ref. [?[ and references 
therein). Several other indirect (often conflicting) hints 
pointing towards the existence of particle DM have been 
proposed over the last few years, and it is therefore im- 
portant to search for strategies that allow to robustly 
and conclusively disentangle a DM signal from a more 
mundane standard astrophysical origin. 

The origin of the cosmic gamma-ray background 
(CGB) measured with EGRET [81] is currently uncertain 
and the the most favored explanation calls for the exis- 
tence of an unresolved population of active galactic nuclei 
(AGNs). Recent determinations of the gamma-ray lu- 
minosity functions (GLF) show however that unresolved 



blazars alone can explain only 20-50% of the measured 
CGB |9j, therefore leaving room for other gamma-ray 
emitters. Besides other standard astrophysical sources, 
e.g. unresolved gamma-ray emission from clusters of 
galaxies [l^, [ll| or normal galaxies [l^, cosmological 
WIMPs annihilation could also contribute to the CGB 

[3,51 EH- 

Assuming a smooth profile for DM halos, the absence 
of intense gamma-ray emission from the center of our 
galaxy constrain the DM contribution from cosmological 
halos to be rather low but it has been shown that 
the presence of substructures can largely boost this signal 
without being in conflict with galactic bounds [l^, fl9| . 
Here we focus on mini-spikes, i.e. large DM overden- 
sities that might form around intermediate-mass black 
holes (IMBHs), due to the adiabatic contraction of the 
DM density profile during the IMBHs' formation and 
growth Unlike the case of a DM spike around the 
central supermassive black hole (SMBH) of our galaxy, 
which would inevitably be disrupte d by the combination 
of several astrophysical processes [l3l [2ll . [22| . the de- 
pletion of mini-spikes is expected to be far less efficient. 
It has been shown that taking into account the contribu- 
tion from cosmological mini-spikes, DM annihilations can 
largely contribute to the measured CGB [l^, while spikes 
around SMBHs can provide only moderate boosts |23|. 

In addition, mini-spikes in the Milky May or nearby 
galaxies such as M31 could be detected with neutrino 
telescopes [2^1, boost anti- matter fluxes 25 1 or individu- 
ally resolved by gamma-ray telescopes such as the Fermi 
LAT [23, |2^ |23| . The simultaneous detection of several 
sources with the same energy spectra, showing a cut-off 
at the DM mass, would be a smoking gun for WIMPs 
annihilation. On the contrary, it is difficult to extract 
straightforward evidences for DM annihilation from the 
study of the CGB spectrum itself. A search for these 
objects based on a HESS survey of the Galactic plane 
region has alread allowed to set some interesting con- 
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straints on the mini-spikes scenario [28[. However, ad- 
ditional information can be extracted by the anisotropy 
data Hi, [33, m m m, [33. in particular, the CGB 
angular power spectrum from blazar and from DM an- 
nihilation in halos or subhalos are quite different, due 
to their different energy spectra, cosmological distribu- 
tion, and radial emissivity profiles. Therefore, the study 
of the CGB angular power spectrum provides, in princi- 
ple, a robust and direct tool to discriminate between the 
two different scenarios. Assuming the unresolved blazar 
contribution as a "known" background, DM annihilation 
could be detected with roughly 2 years of Fermi data, 
provided they contribute a fraction > 0.3 of the CGB at 
10 GeV 

In this paper we perform the angular anisotropy anal- 
ysis for the case of cosmological DM mini-spikes around 
black holes. We compute the angular power spectrum 
for different DM benchmark setups, varying the parti- 
cle mass and the annihilation channel, and for different 
gamma-ray energies, showing that the results are quite 
sensitive to all of these variables. We also discuss the 
possibility to distinguish with Fermi data the mini-spike 
scenario from the case of substructure-dominated emis- 
sion. 

The remainder of the paper is organized as follows: 
Section |TT] is devoted to a discussion of IMBH forma- 
tion and to a summary of the results of recent numer- 
ical simulations. In Sec. IIIIl we compute the contribu- 
tion to the CGB mean intensity from DM annihilation 
in cosmological mini-spikes and from blazars, while in 
Sec. IIVI we compute the angular power spectrum for the 
two cases. A mixed scenario is presented in Sec.|Vl where 
we also discuss prospects for detecting DM annihilation 
with the Fermi Telescope and the effect of changing parti- 
cle DM parameters and the gamma-ray energy at which 
the anisotropy is studied. Finally, conclusions are pre- 
sented in Sec. IVII Throughout this paper, we adopt a 
flat ACDM model with the cosmological parameters from 
WMAP 5-year data [H. 



II. INTERMEDIATE MASS BLACK HOLES 

A. IMBH formation 

Black holes in the range 20 < MbhjMQ < 10^, are 
commonly dubbed as IMBHs (see Ref. [3^] for a review) . 
Despite the lack of conclusive observational evidence for 
the existence of IMBH, many clues have been accumu- 
lated during the last few years, among which the most 
significant is perhaps the detection of ultra luminous X- 
ray sources (ULXs), interpreted in terms of accreting 
IMBHs [13, [111 . Although alternative explanations have 
also been proposed (39l.l4(]|. they seem to be problem- 
atic or in some cases ruled out [H, [4l| , suggesting that 
at least a fraction of the observed ULX sources are in- 
deed IMBHs. Furthermore, the presence of SMBHs at 
high redshifts, as indicated by the detection of high red- 



shift quasars at z '-^ 6 in the Sloan Digital Sky Sur- 
vey [43 . l43l . [4^ , favors a formation scenario with rapid 
accretion and mergers from massive seed BHs [45| . Theo- 
retically, this hierarchical scheme would also help explain 
the tight observed correlations between the SMBH mass 
and the host galaxy properties. IMBHs themselves could 
be the remnants of first stars, also referred to as Pop 
HI stars [s^. These stars, the first formed in the Uni- 
verse, are generally predicted to be very massive, as a 
consequence of their extremely low metallicity and the 
absence of wind and pulsations. While the fate of Pop 
HI stars with masses in the range 100 < M/Mq < 250 is 
perhaps to explode due to pair instability, and to leave no 
compact remnant, heavier stars finally collapse directly 
to IMBHs, with little mass loss. Reference [20*1 investi- 
gated this formation scenario, by following the evolutions 
of the BHs host halos and taking into account IMBH 
mergers, and concluded that at z = a Milky- Way-sized 
halo should host a population of roughly 1000 unmerged 
IMBHs, with masses in the range 10^ - 10^ M©. 

Here we consider an alternative IMBH formation sce- 
nario, pro pos ed in Ref. [46j and further investigated in 
Refs. [20I, [43]. In this framework, IMBHs are formed 
from gas collapsing in DM mini-halos at high redshifts. 
In massive enough halos, the molecular hydrogen cooling 
is very efficient and a proto-galactic disk is formed at the 
center. Gravitational instabilities in the disk then trigger 
an effective viscosity that drives an inflow of the gas lying 
in the low tail of the angular momentum distribution. In 
the absence of halo mergers, the process continues until 
the explosion of the first generation of supernovae, which 
heats up the disk. The mass transferred to the center 
of the halo undergoes gravitational collapse and a BH is 
then rapidly formed. The condition that the BHs forma- 
tion timescale be shorter than the major merger timescale 
and that enough molecular hydrogen be present to form 
a pressure support disk, sets a lower limit on the mass 
of the host halo, i.e., Mcr — 10® Mq at the reionization 
redshift Zre- The IMBH mass function is predicted to be 
a log-normal, with a peak at Mbh — 2.1 x 10^ AIq and 
a spread <7bh = 0.9. 

Reference [l^j studied the population of IMBHs that 
would result from the aforementioned formation scenario 
in our own Galaxy. Specifically, the authors simulated 
the formation of a Milky- Way-like DM halo starting from 
mini-halos at high redshifts, following the hierarchical 
merger history of the latter until z = in the context of 
a ACDM model for structure formation. In that analysis, 
the formation of IMBHs in a given halo follows the pre- 
scription given in Ref. [46] , and pair BH mergers occur if 
the pair distance is lower than 1 kpc. 

IMBH formation is absent after reionization, z < Zre, 
since most of the molecular hydrogen, the main baryonic 
coolant, is ionized. In the simulation of Ref. [13], the 
authors find that IMBHs formation is highly suppressed 
for z > Zre since the suitable hosts for BH formation 
become increasin gly r are as redshift increases. Therefore, 
according to Ref. [43] , the formation redshift distribution 
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Figure 1: Radial distribution of the IMBH population in the 
Milky Way from the numerical results of Ref. [20]. The points 
refer to an average among the 200 Monte Carlo realizations 
of the Milky Way halo and error bars show the scatter among 
realizations. The solid line is an analytical fit and the dotted 
line is a NFW profile. 



is peaked at Zre- 

The authors performed 200 statistical realizations of 
the IMBH population providing for each IMBH its dis- 
tance from the center of the galaxy, its mass and the 
surrounding DM distributions. The average number of 
unmerged IMBHs is Nbh = 101 ± 22. The radial dis- 
tribution of the IMBH population is described by the 
volume probability, g{r), shown in Fig. [T] for an average 
realization among the 200 realizations of the IMBH pop- 
ulation in the Milky Way. The function g(r) is simply 
defined as the probability to find an IMBH at a radial 
distance r from the Galactic center, in a spherical shell 
of thickness dr. The volume probability function is nor- 
malized to 1 between 1 kpc and the maximal distance 
from the Galactic center at which an IMBH is found, 
i.e. roughly 300 kpc. The error bars in the plot reflect the 
scatter among the 200 realizations. 

The distribution is well fitted by the analytical function 



g{r) = 5.9610" 



9.1 kpc 



0.51 



-10.8 



kpc 



The logarithmic slope, 7 = dlog (7/dlogr, is 1.5 at 1 
kpc and 4.5 at 200 kpc, and therefore the resulting dis- 
tribution is cuspier than a Navarro-Frenk- White profile 
(NFW) [11], shown in Fig. [1] for comparison. 

Starting form the simulation of the IMBH population 
in the Milky Way halo in Ref. [20|, the catalog can be 
adapted for other galaxies by rescaling the total number 



of objects by the ratio between the host halo masses, and 
the galactocentric distance by the ratio of virial radii, 
K. The volume probability distribution function of the 
IMBH population in a given galaxy, w(r), is therefore ob- 
tained from that for the Milky Way, as u{r) = g{r/K)/K^. 
This procedure has been satisfactory tested by compar- 
ing the results with a limited number of mock catalogs, 
obtained as an exploratory study in Ref. ,20], for differ- 
ent host galaxies masses. The virial radius for a halo 
of a given mass M at a given redshift is defined as the 
radius of a spherical volume within which the mean den- 
sity is Ac{z) times the critical density at that redshift, 
M — 47rrJ5j^Ac(z)/9c(^3, with the virial overdensity 
Ac(z) as given in Ref. |49l| . 



B. DM distribution around IMBHs 

The process of adiabatic growth of a black hole at the 
center of a DM mini-halo produces a steepening of the ini- 
tial DM profile, leading to large DM overdensities called 
spikes i50;] . Interestingly, while for the SMBHs at the cen- 
ter of the galaxies the DM interaction with baryons and 
dynamical processes tend to weaken the spike [13, l2^ . 
these effects are not effective for the dense mini-spikes 
around IMBHs. Starting from a DM profile with a power 
law behavior in the proximity of the BH p ~ r"'' , the fi- 
nal spike profile after adiabatic growth reads [l^l : 



,(r) = p{rsp) 



sp 



-Isp 



(1) 



where p is the density function of the initial profile and 



the final slope "fsp is given by 7sp 



_ 9-27 
4—7 



, weakly depend- 



ing on its initial value 7. The radius of the gravitational 
infiuence of the black hole rn sets a limit within which 
the spike profile is valid, i.e. r^p « 0.2r;i, where is 
implicitly defined as: 



M{r<rh)= I p{r)r^dr = 2M,, 
Jo 

with M, is the mass of the black hole |22] . 

In the simulation of Ref. [2^ the authors consid- 
ered an initial NFW DM density profile, with the av- 
erage parameters for the spike set as rgp = 6.8 pc and 
Psp — 1.2 10^° Mq kpc^"^. We employ here these refer- 
ence values throughout our analysis. 

After the formation of the BH, the DM number den- 
sity decreases becuase of DM pair annihilations as: — 
— (cru)n^ with (av) the annihilation cross section times 
velocity. The solution to this equation gives an upper 
limit to the DM density pHm = x {(jv)~^{t — tf)~^ 
where indicates the DM mass and t — tf is the time 
elapsed since BH formation. We denote rum the radius 
where this maximum density is reached. The density 
is considered to be constant within a cut-radius defined 
as Tcut = Max[4i?5c/iu), J";™] where Rschw is the BH 
Schwarzschild radius. For the spikes in the simulation 
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Figure 2: Photon spectra for DM annihilation into bb and 
T~^T~ . The DM particle mass is set to = 100 GeV 

of Ref. [13], the cut-radius at z = is, averaging over 
all mini-spikes in the simulation, Vcut = 5 x 10"'^ pc, for 
(av) = 3 10-26 cm3 s'^ 

III. COSMIC GAMMA-RAY BACKGROUND 

We compute in this Section the mean intensity of the 
CCB from unresolved mini-spikes and from unresolved 
blazars. In this and following sections, we use a similar 
notation to that adopted in Ref. [30l |. 

A. Dark Matter Annihilations 

1. DM annihilation contribution to cosmic gamma-ray 
background 

Following Ref. [3] > the CGB gamma-ray flux from cos- 
mological DM mini-spikes, defined as the number of pho- 
tons per unit area, time, solid angle and energy, is then 
obtained as: 

{I{E)dm) ^ j drW {E[l + z], z) , (2) 

where 

W{E, z) = {E[l + z]) e-(-[^+^I-)A2(z). 

(3) 

The absorbtion of gamma-rays due to interaction with 
the diffuse extragalactic background light is parametrized 



through the effective optical depth r as in Ref. [3] • The 
comoving distance r and the redshift z are interchange- 
ably used and the element dr is simply dr = c/H{z)dz 
with H{z) the Hubble function. The function A'^{z) in 
Eq. dSl) is 

A2(z) = n(z) r" p%{r)d\ 

J Tcut 

with n(z) the comoving number density of IMBHs. 

The gamma-ray annihilation spectrum dN~JdE de- 
pends on the DM particle physics model; i.e., it deter- 
mines the branching ratios for annihilation in Standard 
Model final states. Given these branching ratios (which 
can be computed for any specified particle DM model), 
the quantity dN^/dE can be reconstructed via Monte 
Carlo simulations. This is how, for instance, dN^/dE is 
computed in codes like DarkSUSY ^1] which, in partic- 
ular, makes use of Pythia [s^, Hll Monte Carlo simula- 
tions. 

From the discussion above, it is clear that the spe- 
cific DM annihilation spectrum depends critically on the 
particle physics model. In the present study we wish to 
consider a particle dark matter setup as model indepen- 
dent as possible. As such, we consider two representa- 
tive standard model final states, and assume that the 
DM particle annihilates 100% of the time in one of those 
two final states. For definiteness, we consider the final 
states bb and r'*"r~. The choice is motivated by both 
theoretical and phenomenological considerations: first, 
in the context of supersymmetry, perhaps the best mo- 
tivated extension to the standard model encompassing a 
DM candidate, these final states are ubiquitous; second, 
the resulting DM annihilation spectra dN^/dE cover the 
two extreme cases of a soft photon spectrum (bb) and of 
a relatively hard spectrum {t'^t~). Even harder photon 
spectra are in principle possible, for instance in the con- 
text of universal extra dimensions [3| , or in other models 
with a large branching ratio in charged leptons. This is 
not critical to us, since we only focus on a single gamma- 
ray energy in our analysis; our results for the r+r^ are 
conservative with respect to even harder photon spectra, 
and the comparison with the soft spectrum we picked is a 
solid guideline to what would change with an even harder 
spectrum. 

In supersymmetry, in the large tan f3 regime favored by 
Higgs searches at LEP, the dominant annihilation final 
states for the lightest neutralino include gauge bosons (if 
kinematically open) and down-type fermion-antifermion 
final states. The role of gauge bosons depends on the hig- 
gsino fraction of the lightest neutralino. Supersymmet- 
ric models with radiative electroweak symmetry breaking 
and gaugino unification at the grand unification scale fea- 
ture generically a small higgsino fraction. In any case, 
the spectrum resulting from gauge boson final states re- 
sembles closely the bb spectrum d, [H, HI] . If down-type 
fermion-antifermion final states dominate, pair annihila- 
tion into bb is the dominant channel, possibly competing 
with T~^T~ but winning over it by a factor 3 from color 
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Figure 3: Cosmic gamma-ray background spectrum from DM 
annihilation in mini-spikes (dashed) and blazars with the best 
fit LDDE GLF model(dotted). The sum of the two signals is 
shown as a solid line and the data points are from EGRET 
data [Hi. 



and by the square of the bottom-to-tau mass ratio (see 
e.g. [56|). In some cases, however, supersymmetry pre- 
dicts a large branching ratio in t^t^ , for instance when 
the lightest neutralino relic abundance is driven by coan- 
nihilation with the lightest stau, which then also medi- 
ates the dominant pair-annihilation channel. Several su- 
persymmetric models feature r+r" as the dominant an- 
nihilation channel. In addition, other models [3| where 
for instance the quantum numbers of the DM particle 
weigh favorably charged leptons over quarks, naturally 
feature a hard photon spectrum, close to T+r~. 

In summary, in the present study we restrict our- 
selves to the two final states bb and r+r" as represen- 
tative WIMP annihilation final states bracketing a wide 
range of model-dependent predictions. The input spec- 
tra, shown in Fig.O are the results of the numerical study 
of Ref. [53. 



2. Modeling the mini-spikes number density 

The IMBH number densit y is parametrized following 
the numerical study of Ref. [43|. We consider that one 
IMBH is formed at redshift Zre in every DM halo with 
a mass higher than M,nin — IO^Mq. As mentioned in 
Sec. ini IMBHs cannot be formed at more recent epochs 
and the formation at higher redshifts is negligible. In 
any case, the prescription we have adopted provides a 
lower limit to the CGB mean flux from mini-spikes, also 



in view of the fact that more than one IMBH could be 
formed in larger halos as well. 

The comoving number density at the formation red- 
shift, is obtained as: 



n{Zre) 



(4) 



We employ here the halo mass function dn/dM(M, z) 
given in Ref. [s^ , with the transfer function of Ref. [13] • 
After formation, IMBHs get redistributed in halos dur- 
ing their hierarchical mergers. At the present epoch, the 
comoving number density of unmerged IMBHs is given 
by: 



n(0) 



M„ 



dM^iM,^-0)N,, ^^^J^_^^^ , (5) 



with the average number of IMBHs in the Milky Way halo 
Nbh obtained from the simulation of Ref. [43| . Here we 
assume a linear dependence of the number of unmerged 
IMBHs on their host halo mass. As noticed in Ref. 
reasonable deviations from the this linear interpolation 
produce small changes on the final CGB flux. 

At intermediate redshift, we follow the prescription of 
Ref. p^ . and compute n{z) assuming a redshift power- 
law behavior, with the index (3 obtained by fitting n{z) 
at z = and z — Zf : 



n{z) = n{zf) 



1 + z 

l + Zf 



/3 



(6) 



Reference [43] found that a Milky- Way like galaxy would 
host a population of Ngp = 101 IMBHs at z = 0. For 
the same choice of the cosmological parameters and us- 
ing Eqs. HU and (O we obtain (3 = 0.3, as in Ref. [ll|. 
This computation can be updated by using the more pre- 
cise measurements of the cosmological parameters from 
WMAP5 fsE]. Keeping /3 = 0.3, we obtain a sensible de- 
crease of the IMBH number density and for a Milky- Way 
like halo at z—Q, we get Ngp — 40. 

For the rest of the paper we will therefore assume 
Nsp = 40 and /? = 0.3 to parametrize the IMBH number 
density. Using Eq. ([2]) we can now compute the mean 
extragalactic gamma-ray flux from DM annihilation in 
cosmological mini-spikes. The integration over z is per- 
formed up to the formation redshift, i.e., z^e- The re- 
sults are shown in Fig. [3] adopting = 100 GeV and 
{av) = 3 X lO^'^^ cm^'^ s~^ and for DM annihilation into 
bb and t'^t^ . In the same plot are shown the measure- 
ments of the CGB extracted from EGRET data 58]. The 
predictions largely depend on the annihilation spectrum, 
with the CGB flux peaking at higher energies for harder 
spectra. For energies of the order 0(1-10) GeV, the con- 
tribution from DM annihilation is at the same level of 
the CGB intensity inferred from EGRET measurement, 
suggesting therefore that in this energy range DM anni- 
hilation could substantially contribute to the total CGB 
flux. 
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In our analysis, we have included also the contri- 
bution from low redshifts, where IMBHs are poten- 
tially detectable. Previous studies have shown that the 
Fermi satellite is expected to resolve mini-spikes in our 
galaxy [20|| and maybe Andromeda [2f?| but not further. 
On the other hand, the contribution of IMBHs from 
z < 10~^ to the extragalactic gamma-ray background 
is negligible. 



B. Unresolved Blazars 

The gamma-ray luminosity function (GLF) of blazars 
is obtained from the luminosity dependent density evo- 
lution (LDDE) model of Ref. |6i|. 

The CGB flux from unresolved blazar is computed as: 

E{Ib{E))=J dzj-^J^ dLp^{L,z)TE{L,z). 

The functions in the Equation above are derived in 
Ref. fsoj and references therein. The minimum blazar 
luminosity is taken to be L,„i„ — 10^^ erg s^^ and the 
EGRET flux sensitivity above 100 MeV is lO""^ cm^^ 
s~^. In Fig. [3] we show the results for the best-fit LDDE 
GLF model (details on the blazar model can be found in 
Ref. 01). 

IV. COSMIC GAMMA-RAY ANGULAR 
CORRELATIONS 

A. Dark Matter Annihilations 

The angular power spectrum Ci of the CGB from 
DM annihilation in substructures has been computed in 
Ref. [30]. Here we adapt their formalism to the case of 
mini-spikes and we refer to the original reference for the 
derivation of the equations. 

The angular power spectrum from mini-spikes is ob- 
tained as: 

{I{E)fCi =J^W ([1 + z]E, zfPoM {^-.z^ , (7) 

where PoMik) is the spatial power spectrum of mini- 
spikes and it can be divided into 1-halo and 2-halo terms: 

PDM{k) - P^\k) + P^\k), (8) 

^'"<^'^£/"^(^)'|»<*'"'|"' 

X P''"'='"'(fc,z). (10) 

These terms refer to correlations between two points in 
the same halo (1-halo) or in two different halos (2-halo). 
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Figure 4: Angular power spectrum from DM annihilation 
in mini-spikes computed at i?o = 10 GeV. We show sepa- 
rately the contributions from 1-halo (red) and 2-halo (blue) 
terms. The total angular power spectrum is the sum of the 
two curves. Solid lines refer to DM annihilation into 66 and 
dashed ones are for t^t~ final states. We set the DM particle 
mass to = 100 GeV. 

The function w(fc, M) is the Fourier transform of the 
IMBH volume probability, defined in Sec. |TT1 The lin- 
ear power spectrum P(fc) is obtained using the transfer 
function of Ref. [g^I and the bias parameter is taken from 
Ref. 62]. 

The function {N\M) gives the number of IMBHs in an 
halo of given mass at a given redshift and it is related to 
the IMBH comoving number density as: 

n(z)=^J dM^{N\M}. (11) 

As noticed in Sec. |TT1 at z = {N\M) is well approxi- 
mated by {N\M)un = A^sp doia-i^-iMg ) ^i*'^ ^^v corre- 
sponding to A^sp — 40, as appropriate for the Milky Way. 
On the other hand, at the formation redshift we assume 
that one BH is formed for every halo with mass above 
Mmin - Formally, there is no unique expression for {N\M) 
which interpolates the two regimes above and that, at the 
same time, allows one to reproduce Eq. ([5]) from Eq. (jlip . 
We have explored different parametrization for (A^|Af) 
encompassing its limiting behaviors at z = and z = Zre 
and overestimating and underestimating n(z) with re- 
spect to Eq. ([5]). Since the power spectrum computed 
in Eq. ([7]) is dominated by the contribution at small z, 
we have found that these different choices produce dif- 
ferences in Ci always between a factor 2, that are within 
other uncertainties in the calculations. This is also true 
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also for the cross-correlation terms that we will introduce 
in Sec. |Vl 

From Eq. ([7]), we note that the multipoles C/ are in- 
dependent of the value of (cru) and of the choice of DM 
density profile around each IMBH. We also find that they 
are weakly dependent on the normalization iVb/j. 

In Fig. |4] we show, for the two different WIMP an- 
nihilation channels, the contributions of 1-halo and 2- 
halo terms on the angular power spectrum. We picked a 
gamma-ray energy at which we compute the anisotropy 
power spectrum of i? = 10 GeV, and fixed the particle 
DM mass to = 100 GeV. 

The 2-halo term turns out to be negligible at all an- 
gular scales. The slope of the 1-halo term lies between 
those of the 1-halo terms for annihilation in subhalos 
and smooth NFW halos computed in Ref. [30] (see their 
Fig. 2). This can be understood considering that the 
signal in the subhalo- and in the smooth-halo-dominated 
cases follow respectively the density profile and its square 
and for the case of a NFW profile the two distributions 
are respectively steeper and shallower than the IMBH 
radial distribution. The increased normalization of the 
power spectrum with respect to the case of subhalos emis- 
sion is explained by the same argument: the Fourier 
transform of the IMBH profile gets more power at high k 
with respect to that of NFW. The same tendency is found 
for the two choices of the annihilation spectra. The angu- 
lar power spectrum for DM annihilation into bb is larger 
than that for r+r^ final states because at the energy of 
Eq = 10 GeV, the former photon spectrum is significantly 
steeper than the latter. 



B. Blazars 



The angular power spectrum from unresolved blazar 
comes from the contributions of a Poisson term Cf and 
a correlation term Cp, respectively the 1-halo and 2-halo 
terms: 
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The blazar bias 6b indicates how strong blazars are clus- 
tered with compared to the linear matter power spec- 
trum. Presently, this value is uncertain, and different 
(generically inconsistent) estimates are inferred from dif- 
ferent techniques. Current observations give an upper 
bound 6s < 5 (see Ref. [61j]). 
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Figure 5: Angular power spectrum of the CGB from unre- 
solved blazars expected for Fermi. We show separately the 
Poisson (dotted) and the correlation (dashed) terms. The to- 
tal is simply the sum, and is shown as a solid curve. We 
assume here the best-fit LDDE GLF model. 



Reference [30| estimated the correlation term assuming 
either a bias model inferred from quasar observations or 
a simply constant 6b = 1. The results obtained are quite 
similar since the main contribution to the CGB comes 
from low-redshift blazars, which have bias close to 1. In 
addition, for I > 10 the total angular power spectrum is 
dominated by the Poisson term. 

We present in Fig. [5] our predictions for the angular 
power spectrum expected to be reconstructed from Fermi 
data, adopting the best-fit LDDE GLF model. We as- 
sume a Fermi point source sensitivity of 2 x 10~^ cm~^ 
s~^, the value expected for energies above E = 100 MeV 
and two years of full sky survey mode, for sources with 
a spectral index equals to 2. We note that the power 
spectrum is independent of the gamma-ray energy, since 
we have assumed the same power-law spectrum for all 
blazars and these dependence exactly cancels when we 
divide by the mean intensity squared in Eq. (fT3|) and 
Eq. nil). 



V. DISTINGUISHING DARK MATTER 
ANNIHILATION FROM BLAZARS 



We outline here the prospects for distinguishing DM 
annihilation from blazar emissions in the angular power 
spectrum of CGB with Fermi. 
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Figure 6: Angular power spectrum of the CGB from DM annihilations around IMBHs at a photon energy Eq = 10 GeV. 
Dashed line shows the contribution from DM annihilation {foMCP'^'), dotted line is for blazars {/bCP) and the dot-dashed 
line is the cross-correlation term 2fDMfBCi'^. The total signal Cf is shown as a thick black solid curve. Error bars are for 
2-years of Fermi data. The thin blue solid curve show the DM signal for DM annihilations in subhalos [sO] (see text for more 
details) . 



A. Angular correlations of CGB in the two 
component case 

The CGB background receives contributions from 
both DM annihilation and from ordinary astrophysical 
sources, with unresolved blazars being a representative 
candidate for the latter class of emitters. For the detec- 
tion of DM annihilation in the CGB, blazars therefore 
constitute a background. Their contribution is currently 
uncertain but we expect it will be modeled rather pre- 
cisely with the Fermi catalog of detected blazars. In ad- 
dition, as mentioned in Sec. llVl the angular power spec- 
trum for astrophysical sources is energy independent and 
therefore it could be calibrated at low energies where 
the contribution from DM annihilation is negligible and 
then subtracted from the total anisotropy data. For this 
analysis we therefore treat the blazar contribution as a 
known background, and we study the prospects for de- 
tecting DM annihilation on top of it. 

In this two component analysis, the total CGB inten- 
sity is the sum of the DM and blazar contributions: 

{IcGBm = {Idm{E)) + {Ib{E)). 

Labeling with /dm the fraction of the total CGB coming 
from DM annihilation, foM = {Idm{E)) / {Icgb{E)) , 
the total angular power spectrum is: 



_ f2 

— JDM'^l 



2fDMfBCf'' + fgCi 



p2 r-iB 



where Cj and C; are respectively the angular power 



spectrum from DM annihilation and blazars and Jb is 
simply /b = 1- foM- 

The cross-correlation term Cf^ has been studied in 
Ref. [301 and is divided into 1-halo and 2-halo terms: 



a 



Cr.l — halo 



W{[l + z]E,z) 
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A relation between the blazar luminosity and its host 
halo mass, M[L\ is given in Ref. 61]. 

In this two component framework, the total signal Cf 
and the background noise Cf therefore read: 



= foM + 2/dm(1 - fDM)Ci 
= (1-/. 



DM) '^l 



(17) 
(18) 



The GLF-LDDE blazar model in Ref. ^ basically 
depends on three parameters (71,(7, fc) and as reminded 
in Sec. Illli the best-fit model only accounts for ~15% 
of the CGB intensity at 10 GeV. However, varying the 
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parameters of the blazar model allows to explain differ- 
ent fractions of the CGB. For example setting them to 
(71 = 1.36, q = 3.80, k = 3.15 x 10"*^) we obtain a blazar 
fraction Jb — 0.71. On the other hand, the contribution 
from DM annihilation in mini-spikes is largely affected 
by astrophysical and particle physics uncertainties. For 
example, in WIMP models the mass usually lies in the 
broad range 0(1-1000) GeV^ and {<tv) can largely differ 
from the thermal value {av) = 3 x 10~^^ cm^^ s~^ in the 
presence of efficient coannihilations or Sommerfeld cor- 
rections, if the DM candidate is nonthermally produced 
or if a modified cosmological expansion rate is postulated 
at the time of WIMP freeze-out. Moreover, as discussed 
in Sec. IIIH the number of mini-spikes in halos could dif- 
fer from that found in simulations, since IMBH forma- 
tion could have been underestimated or on the contrary 
IMBHs could have been more efficiently destroyed by as- 
trophysical processes than what is expected. In addition, 
the DM density profile around each IMBH could be mod- 
ified as well by feedback. 

Motivated by these arguments we compute two differ- 
ent models of blazars, explaining respectively a small and 
an high fraction of the CGB at 10 GeV, and we assume 
that the remaining CGB intensity comes from DM an- 
nihilation in mini-spikes. As a benchmark model we fix 
the DM mass to — 100 GeV and we refer to anni- 
hilation to bb pairs. Following Ref. [sOl, we choose an 
energy of observation Eq = 10 GeV as a compromise 
between maximization of signal count and minimization 
of the Galactic emission. At lower energies, the galactic 
foreground becomes stronger, masquerading the extra- 
galactic component, while at higher energies, the pho- 
ton count is more suppressed. However, we perform our 
analysis also for different choices of DM parameters and 
energies of detection. 



B. Prospect for detection with the Fermi Telescope 

The Large Area Telescope (LAT) onboard the Fermi 
satellite is currently taking scientific data in a survey 
mode. The LAT has a more than one order of magni- 
tude better sensitivity in the 20 MeV to 10 GeV region 
than its predecessor Energetic Gamma Ray Experimen- 
tal Telescope (EGRET) onboard the Compton Gamma- 
ray Observatory In addition, the LAT extends the 
high-energy gamma-ray region up to around 300 GeV. 
In the present study, we consider a mean exposure of 
1.2 X 10^^ cm^ s, corresponding, roughly, to 2 years of 
all-sky survey mode operation [63, ISa, [63| . We assume 
an angular resolution for 68% containment of the point 
spread function of ai, ~ 0.115°, appropriate for energies 
of around 10 GeV. Our choices reflect those described 



^ For a recent discussion of ultra-light WIMPs in supersymmetry 
see Ref. |63|| . 



in Ref. [30|. The angular resolution improves at larger 
energies, and degrades at lower energies. 

For the type of study hereby presented, a thorough 
knowledge of the gamma-ray galactic background will be 
warranted. In addition, disentangling the diffuse extra- 
galactic background from the mentioned galactic emis- 
sion will also be challenging. Realistically, the 2 years 
of observations we consider refer not to the early stages 
of the mission but, rather, to a stage when these back- 
grounds are considered to be thoroughly under control. 

Considering the Fermi specifications described above, 
the projected 1-a error bars of the CGB power spectrum 
from DM annihilation is: 

We take a bin width Al = 0.51. The window function of 
a gaussian point spread function is Wi = exp{—Pa^/2). 
Cat is the photon spectrum of the photon noise and it 

is given by Cn = ^skyNtot/N^cB '^ith ^tot and Ncgb 
respectively the total and CGB photon numbers detected 
from a region of s ky Qsky ■ 

Following Ref. [30| , we restrict the analysis to galactic 
latitudes \b\ > 20°. At lower latitudes, the galactic fore- 
ground dominates over the CGB flux, while the situation 
is expected to be reversed in the region we consider. Af- 
ter the cut of the galactic plane, the fraction of sky we 
consider is fsky = 0.66. Using Ntot Ncgb we obtain 
Cn ~ 4:nfsky/NcGB = 8 x 10-^{E/10 GeV). Here we 
employ the total CGB flux as estimated from EGRET 
data. We note, however, that since Fermi is expected to 
detect a large number of blazars, the total GCB inten- 
sity will be in all likelihood reduced, possibly lowering 
our error estimations. 

In Fig. [51 we present our predictions for two blazar 
models contributing a fraction /b = 0.13 and fs ~ 0.71 
of the total CGB flux at Eq = 10 GeV. We show the 
signal and the background power spectra that Fermi is 
expected to measure after two years of observations as 
well as the projected 1-cr signal error bars. The signal is 
detected if Cf > (5Cf . We notice that this occurs even if 
the DM contribution is very small. In addition, the shape 
of the DM power spectrum is very different from the one 
corresponding to blazars. This feature could therefore 
help distin guis h the two scenarios. 

In Refs. [29l. [30|. the authors first studied the angular 
anisotropies of the CGB from DM annihilation. They fo- 
cused their attention on two scenarios, the first assuming 
that the DM signal is dominated by annihilations occur- 
ring in cosmological DM halos and the latter considering 
that the dominant contribution comes from the popula- 
tions of DM clumps hosted in the main DM halos. For 
each possibility they took into account the possible un- 
certainties on the minimum halo mass value and on the 
halo occupation distribution, i.e., the number of subha- 
los in a parent halo of given mass. For these frameworks, 
they computed the angular power spectrum from DM 
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annihilation that Fermi is expected to measure. They 
concluded that provided DM annihilation contribute to 
the CGB at 10 GeV with a fraction /dj\/ 0.3, after two 
years of data taking Fermi will be able to detect the DM 
signal. 

In Fig. \6\ we show their results for the most promis- 
ing case, i.e., when the DM signal is dominated by cos- 
mological clumps with an halo occupation distribution 
{N\AI) (X M. We consider that DM annihilations in sub- 
halos contribute to a certain fraction fuM to the CGB 
intensity at 10 GeV and the remaining flux comes from 
blazars. In each plot the DM fraction Jdm is the same 
for the mini-spike and clump scenarios. The signal for 
the subhalo case with the associated 1-cr error bars is 
plotted as a thin blue solid line. Comparing the DM sig- 
nals in the mini-spike and subhalos scenarios we notice 
that provided that DM annihilations largely contribute 
to the CGB mean intensity, there are promising prospects 
for distinguish the two cases. This conclusion is further 
strengthened if we consider DM annihilation in cosmo- 
logical smooth halos instead of clumps, since, as stated 
before, the expected angular power spectrum is smaller 
than when subhalos emission dominates. 



C. Power spectrum dependence on energy of 
detection, annihilation spectrum and DM mass. 

Even if the results discussed in the section above refer 
to a certain specific choice of DM parameters and energy 
of detection, we have also repeated the calculations for 
different cases. Rather than presenting all the plots, we 
just show in Fig. [7| what we obtained for some bench- 
marks and we try to summarize some general guidelines. 
A more complete analysis, for example dedicated to the 
optimization of the energy of detection as a function of 
the particle mass, is beyond the scope of this study. 

First we show how our predictions change if we pick 
another energy of detection. At energies higher than 
10 GeV the galactic foreground is sensibly suppressed 
but also the photon number from DM annihilation is 
reduced, since the interval of integration in energy is 
shrunk. Therefore, it is not trivial to infer which is the 
effect on the DM angular power spectrum and its error 
bars. We find that even if the CGB mean intensity at 20 
GeV is reduced with respect to its value at 10 GeV, the 
power spectrum is increased. We remind that the power 
spectrum is normalized to the mean flux, as in Eq. ([7]). 
At an energy of 1 GeV the CGB mean intensity is com- 
parable with the galactic foreground therefore in Eq. (fT9|) 
we consider Cn ~ "^flsky/NcGB- For this gamma-ray en- 
ergy, the signal is sensibly reduced and the prospects for 
detection are degraded. 

As pointed out in Sec. lIVI for softer energy spectra, the 
normalization of angular power spectrum is decreased. 
We indeed find this behavior when we compare the re- 
sults obtained for the t'^t~ and for the bb DM annihila- 
tion final states. For the "pessimistic" case of annihila- 



tions into r+T" , assuming an energy of detection of 10 
GeV and = 100 GeV, DM annihilations have to con- 
tribute at least with a fraction /da/ ^ 0.3 to the mean 
CGB in order to be detectable in the CGB angular power 
spectrum with Fermi. 

We finally show the results for = 1000 GcV and 
Eq = 10 GeV. The photon spectra from DM annihilation 
can in good approximation be scaled with the particle 
mass defining the adimensional variable x = E/m-^, as 
in Fig. m Therefore, looking at Eq. ^ and Eq. ([7]), we 
note that there is an approximate scaling which links the 
angular power spectra computed at different energies of 
observations and for particle of different masses. For ex- 
ample, the choices {Eq = 10 GeV, — 1000 GeV) and 
{Eq = 1 GeV, = 100 GeV) actually correspond to 
the same value of dNj/dx, and should thus lead to iden- 
tical angular power spectra. This scaling is broken by 
the dependence of the function r(z) on Eq in Eq. ([3]), 
which fortunately is not important at the energies of in- 
terest, and our qualitative considerations are still roughly 
valid. This can be seen noting that the power spectra in 
Fig. [7] for the two cases above are indeed very similar. 
Note however that, as already stressed, different energies 
of observations significantly affect the projected errors 
bars. 



VI. CONCLUSIONS 

DM annihilation in mini-spikes around IMBHs is a 
promising scenario for indirect DM searches with gamma 
rays. In particular, Fermi is expected to detect a sig- 
nificant fraction of the IMBH population in the Milky 
Way and maybe a few sources in the Andromeda galaxy. 
The remaining cosmological mini-spikes will remain un- 
resolved, but could leave their imprint in the CGB. As 
shown in (isj . for a standard neutralino with a mass 
niy, = 100 GeV and a "thermal" annihilation cross sec- 
tion (ctw) = 3 10~^^ cm^s~^ the predicted CGB flux from 
cosmological mini-spikes is comparable to the EGRET 
CGB flux at gamma-ray energies of 0{1 — 10) GeV. We 
find that, for example, this corresponds to a fraction 
foM = 0.35 and Jdm = 0.72 of the CGB at £; = 10 
GeV, respectively for DM annihilation into bb and r+r^. 
Fermi is expected to resolve a much larger number of 
galactic and extragalactic gamma-ray sources compared 
to its predecessor EGRET, with the expectation of re- 
ducing the measured unresolved diffuse CGB flux. At 
the same time, only IMBHs very close to us will be re- 
solved, therefore the DM contribution to the CGB could 
be increased with respect to our estimates, based on the 
mean CGB spectrum extracted from EGRET data. 

However, in absence of characteristic spectral features, 
it will be problematic to distinguish DM annihilation and 
ordinary astrophysical emissions from the mean GCB 
intensity. Instead, Ref. [30l | showed that gamma-ray 
anisotropy data could provide a more suitable tool to 
pursue this program. In fact, provided DM annihilation 
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Figure 7; Angular power spectrum of the CGB from DM annihilation and for blazars. Lines are as in Fig. |6l Annihilation 
channel, energy of detection, DM mass and fractional contribution to the CGB mean intensity are specified for each panel. 



contributes substantially to the CGB mean intensity, it 
will be detectable in the CGB angular power spectrum 
by Fermi. 

Motivated by these considerations, in this paper, we 
studied the anisotropics of the CGB in the mini-spikes 
scenario. Astrophysical and particle physics uncertain- 
ties largely affect the predictions for the mean CGB in- 
tensity from mini-spikes and also the blazar contribution 
is currently unknown. Considering these two sources as 
the main components of the CGB, we computed their 
angular power spectra for different relative contributions 



and, treating the blazar component as a known back- 
ground, we studied the prospects for DM annihilation 
detection in the CGB angular power spectrum with two 
years of Fermi observations. We expect that consider- 
ing unresolved blazars as a background is a reasonable 
assumption, since their GLF and bias should be quite 
reliably reconstructed from the Fermi source catalog. 

We repeated our computations for different detection 
energies, particle masses and annihilation modes, show- 
ing that our results are significantly affected by all these 
parameters. Interestingly, this could mean that informa- 
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tion on these three quantities can actuaUy be inferred 
from the measured DM-induced gamma-ray anisotropy 
power spectrum. However a more detailed analysis would 
be necessary to fully study the potential of the anisotropy 
technique to reconstruct these parameters. We found 
that the shape of the DM power spectrum is very dif- 
ferent from that of blazars, providing a robust handle to 
disentangle the two signals. Astrophysical sources other 
than blazars could however also contribute to the CGB 
and, if spatially extended, as clusters of galaxies, the 
shape of their angular power spectrum could significantly 
differ from that of blazars, which is dominated at large 
multipoles by the Poisson term. We stress that even in 
this case we could calibrate the astrophysical power spec- 
trum at low energies, where DM annihilations are negli- 
gible and subtract it from the measured total CGB power 
spectrum at the energies of interest. In fact, for sources 
with power-law energy spectra, the gamma-ray angular 
power spectrum is energy independent and this condition 
is common to almost any class of standard astrophysical 



gamma-ray emitter. 

In conclusion, we showed that the prospects for detect- 
ing DM annihilation from cosmological mini-spikes in the 
angular CGB power spectrum with Fermi are promising, 
and that the analysis of the anisotropy power spectrum 
allows not only a discrimination of a DM component 
against astrophysical sources, but also a better under- 
standing of the structures where the DM signal origi- 
nates. 
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